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SECTION  I 

INTRODUCTIUN 

Bleistein  and  Cohen  have  shown  [1,2J  that  a source  distribution, 
giving  rise  to  fields  that  satisfy  the  inhomogeneous  Helmholtz  wave 
equation  (for  acoustics,  electromagnetics,  etc.),  ‘^an  be  separated  into 
radiating  and  non-radiating  portions;  and  that  the  radiating  source 
distribution  can  be  determined  completely  and  uniquely  in  a well-conditioned 
manner  from  knowledge  of  the  radiated  fields  on  any  closed  surface  exterior 
to  the  support  of  the  source  distribution. 

An  alternative  formulation  of  this  separation  is  developed  in  a 
form  suitable  for  deriving  a wave  equation  which  governs  the  radiating 
portion  of  the  source  distribution. 


SKtTION  II 


AN  ALThRN Al  I Vh  KORMIILAT i UN  OF  I'HI  Bl.LISTfalN  t:OHKN  LCjUAl'IONS 


Let  Y and  be  the  normalized  over  6 t'  (O.n)  and 

nui  1 

>()  G (O.ifl)  spherical  harmonic  and  the  normalized  over  r G 10,")  spherical 
Bessel  functions  respectively.  It  thus  follows  that  the  functions 
Y^i^^(e,((>)  jj^(kr)  are  orthonormal  and  complete  in  all  space  [3J,  and  that 
any  function  can  be  expanded  as 


m=-i  1=0 


where  the  coefficients  given  by 


<*■  Jn  n 

C , - I I I Y (b,(J))  j\(kr)  p(e,(|),r)  sin  0 de  dcp  r''  dr  . (2) 

mni  III  mn  ■^1 

•'o-'o  •'o 

If  the  function  p represents  a source  distribution,  giving  rise  to 
fields  III,  which  satisfy  the  inhomogeneous  reduced  Helmholtz  wave  equation 


ly  - -o 


and  the  support  of  these  sources  p is  confined  to  the  interior  of  a sphere 
of  radius  ci,  'hen  the  radiated  fields  ij/  at  the  exterior  of  this  sphere  are 
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ij;(0,4i,r)  - Ik 


L S f f f 
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Y (0'  ,(>' ) j (Kr' ) 
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fi=C 
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,il' .1  ' ' 3'n  •>'  de*  da*  ' *'’  dr'  , r 


X 


which  by  (2)  can  be  written  as 


n 

^(0,^,r)  = ik  V'^r)  . (5) 

11=0  m=-n 

i.e.,  only  the  diagonal  coefficients  C ,of  the  complete  set  of  coefficients 

num 

^mnl^  expansion  (1,2)  give  rise  to  radiated  fields. 

If  the  functions  G are  defined  as 


mnl 


Y 

mil 


(6,  ij)) 


j^(kr) 


(6) 


then  by  the  general  solution  in  spherical  coordinates  of  the  homogeneous 
reduced  Helmholt?  wave  equation. 


G + k2  G 
mnn  mnn 


0 


(7) 


Thus,  by  Green’s  theorem,  (3),  and  (2) 
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dS*(G  Vi;)  - i))  VG  ) = I dv  (G  V^iL  - ili  V^q  \ 
mnn  ^ mnn  | mim  ^ ^ mnii 
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dv  G ;• 
mnn 


= C 

mnn 


(8) 


i.e.,  the  diagonal  coefficients,  which  completely  and  uniquely  determine 
the  radiating  portion  of  the  source  distribution,  can  be  determined  in  a 


■ 


well -conditioned  manner  from  knowledge  of  the  field  >1/  on  any  closed  surface 
s,  exterior  to  the  support  of  the  source  distribution  p. 


SECTION  III 


A WAVE  EQUATION  t-OR  RAUIATING  SOURCE  DISTRIBUTIONS 


The  expansion  (1)  of  an  arbitrary  source  distribution  ()(e,<|),r)  can 
be  broken  into  the  diagonal  and  non-diagonal  parts 


00  1^00 


p ( 6 


‘X)  Z)  X)  ‘=«nl 


n=0  m=-ii  1=0 
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+ C Y (0.<t,)  j (kri 

mnn  mn  -^n 
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19) 


which  can  be  identified  respectively  as  the  non-radiating  (po)  and 
radiating  (p^)  portion  of  the  source  distribution;  i.e.. 


p(0,4>.r)  = po(0,(}>,r)  + p^(0,()),r)  , (10) 


where 
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(11) 


(12) 


By  the  general  solution  in  spherical  coordinates  of  the  homogeneous 
reduced  Helmholtz  wave  equation,  it  follows  from  (12)  that 


b 


An  alternative,  and  much  simpler,  coordinate  inde])endent  derivation 


of  (13)  for  the  far-field  radiating  source  distribution  can  be  accomplished 
in  the  following  fashion. 

llie  three-dimensional  (spatial)  Fourier  transform  of  any  source 
distribution 


d(k,u))  = po<k,u)  1-  P_^(k, id)  . (17) 

However,  by  (16),  p^(k,iD)  is  the  complete  set  of  the  (range  and 
phase  normalized)  radiated  far-field  of  the  source  distribution  d(x,d), 
and  by  (13),  0o(k,tD)  contributes  nothing  to  the  radiated  far-field.  it 
thus  follows  that 


(5o(k,uj)  d^k 


(18) 
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Po(X,(jli)  = 


(2ii) 


J.-ik-i 

Vkj<- 


is  the  non  far* field  radiating  portion  of  the  source  distribution,  and 


k>  (X,Cli}  *“  ~7 

' tXJ  I 


271) 


0 ( k,u))  d ^k 


(19) 


V k=- 


is  the  far-field  radiating  portion  of  the  source  distribution;  where 


p(X,uj)  - po(X,u))  + p (X,ii)) 


(20) 


It  follows  directly  from  (19)  that 


V^p  + p 


U>  (JO 


0 . k = ^ 


(21) 


As  suggested  by  Bleistein  [4],  expansion  of  the  integrand  term 

— i k • X . 

e in  (14),  (18),  and  (19)  in  spherical  harmonics  reveals  that  the 

radiating  and  non-radiating  source  distributions  of  this  alternative 
derivation  are  indeed  the  same  radiating  and  non-radiating  source  distributions 
given  by  (11)  and  (12)  respectively. 


It  can  similarly  be  shown  that  a vector  current  density  distribution 
J,  which  give^  vise  to  electromagnetic  fields  that  satisfy  Maxwell's 
equations,  an  similarly  be  separated  into  radiating  and  non-radiating 
current  density  distributions,  and  that  the  radiating  current  density 
distribution  satisfies  the  vector  wave  equation 


VxVxJ 
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(22) 
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Source  distributions  that  fail  to  satijfy  the  wave  equation  should 
be  carefully  excluded  from  solutions  to  the  inverse  source,  inverse 
scattering,  and  synthesis  problem.  Such  an  exclusion  might  advantageously 
be  achieved  by  viewing  this  wave  equation  as  an  additional  necessary 
condition  (e.g.,  facilitating  regularization  of  ill-conditioned  and/or 

ill-posed  problems)  on  the  solutions  to  the  inverse  source,  inverse 

,i 

j scattering,  and  synthesis  problem. 
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